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In this paper we propose a framework for simulating approximate thermal particle production 
in condensed matter systems. The procedure we describe can be realized by means of a quantum 
quench of a parameter in the model. 

In order to support this claim, we study quadratic fermionic systems in one and two dimensions 
by means of analytical and numerical techniques. In particular, we are able to show that a class of 
observables associated to Unruh-de Witt detectors are very relevant for this type of setup and that 
exhibit approximate thermalization. 



I. INTRODUCTION 

The relative weakness of the gravitational field makes 
difficult any direct test of some of the phenomena 
that should characterize quantum field theory (QFT) in 
curved spacetimes, namely various instances of particle 
creation from the quantum vacuum. For the same rea- 
son, the regime in which quantum gravitational phenom- 
ena deviate significantly from the predictions of Einstein 
gravity (viewed as the effective field theory description of 
a full-fledged quantum gravity theory) seems to be still 
out of reach. Therefore simulators, i.e. models offering 
analogies with gravitational systems, could help us to 
develop concepts, methods and insights for the investiga- 
tion of the relationship between quantum mechanics and 
the geometry of spacetime. 

Since the seminal work of Unruh [1[ , who showed how 
four dimensional Lorentzian geometries can be simulated 
in hydrodynamical systems, the area of analogue mod- 
els has flourished [51, providing nowadays a rather large 
number of (classical and quantum) models based on con- 
densed matter systems, which are able to simulate QFT 
in curved spaces. Furthermore, we have by now numeri- 
cal and experimental results [3HS] which allow us to dis- 
cuss concretely the validity of the theoretical predictions 
for phenomena like Hawking radiation and cosmological 
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particle creation^. 

Most of the analogues are based on continuous approx- 
imations, hence applying methods of continuuum field 
theories. However, this is not the only possibility that 
we have. It has been already proposed to use lattice sys- 
tems to model new interesting situations. For instance, 
the use of optical lattices has been advocated in [TD] and 
[111 I12j in order to explore regimes in which conventional 
continuum mean field descriptions fail or in which the 
quantum vacuum is explored in the strong external field 
limit, regimes difficult to achieve and to control in other 
quantum systems like BECs^. 

Another motivation for discrete analogues comes from 
the number of discrete Quantum Gravity proposals. In 
particular, the possibility that at a microscopic level con- 
tinuum spacetime might be replaced by a completely 
combinatorial structure in the shape of a random graph 
(variously decorated with geometric data) is a widespread 
idea [13]. Despite the variety of proposals, the manipu- 
lation of random graphs (or random complexes) within 
a reparametrization invariant framework presents diffi- 
culties which are so far not completely under control, 
especially if we try to obtain a quantitative explanation 



^ The crucial point is the fact that these phenomena do not need 
Einstein's equations for the metric to be satisfied. They arc in a 
sense kinematical, as far as the metric tensor is concerned. 

^ For instance, the curved acoustic geometry of BECs is related to 
the velocity of the condensate flow, its density, and derivatives. 
To achieve strong gravity efi'ects, then, requires manipulations of 
these quantities without destroying the condensate regime. 
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of the emergence of classical continuum spacetime and 
its dynamics. 

In fact, in |T4HT7] it has been proposed to address a 
simpler problem, that is the emergence of space alone, as 
the outcome of the collective dynamics a certain specific 
class of (quantum) Hamiltonians for (quantum) graphs: 
Quantum Graphity. Since in such a framework there 
is a preferred foliation related to the presence of a pre- 
ferred notion of time (here radically different from space) , 
the inclusion of Lorentz invariance (and of the full four 
dimensional diffeoniorphism invariance) is problematic'^. 
Nonetheless, they represent ideal candidates to discuss 
in an explicit way at least some general features of the 
otherwise vague idea that space, time and geometry are 
emerging from pregeometric degrees of freedom. For such 
a purpose, different models have been introduced so far, 
with the objective of elucidating the way in which the 
dimensionality of the graphs is dynamically controlled, 
as well as the properties of the phase transition from the 
disordered, pregeometric phase to the ordered, geometric 
one [T7| [20l [21] . Possible cosmological effects due to their 
fundamental discrete structures have been considered re- 
cently in m]. 

As we shall argue in detail in the rest of the paper, 
these systems might still offer interesting insights as far 
as quantum field theory in curved spaces is concerned. In 
particular, we will show how these models can be used to 
simulate particle creation phenomena (either in a black 
hole- like configuration, like in [23], or in a cosmological 
setting) by analogy with Unruh effect, tying it to quench 
experiments, which are indeed the natural experimental 
setup. 

Concretely, we will work with a Fermi-Hubbard model 
on a fixed lattice with an additional term in the Hamil- 
tonian, which does not conserve the particle number. A 
similar model appeared already in the literature of Lieb- 
Robinson bounds and has been experimentally probed 

m- 

The model which we consider has been initially derived 
from Quantum Graphity [53], as an effective description 
in which the graph is essentially frozen in a given config- 
uration. The derivation is explained in Section IV in the 
section dedicated to Quantum Graphity. The application 
of our setup to quantum graphity's trapped surfaces, for 
the sake of clarity, is at the very end of the paper, as in 
principle this setup could apply to other physical models. 

The paper is organized as follows. Section II is devoted 
to the detailed explanation of the model and of the calcu- 
lations that will be presented, in particular with respect 



to the notion of particle and detectors. In section III 
we first show the numerical results on the thermality of 
the detected particles, for various 2-dimensional square- 
shaped lattices and for the 1-dimensional ring. The latter 
is treated analytically in full detail. A subsection is de- 
voted to the clarification of the results from the point of 
view of quantum quenches. Section IV explains the ap- 
plication to quantum graphity's trapped surfaces, with 
emphasis on the subtle but important differences with 
the case of true black holes. Some final remarks conclude 
the paper. 

II. THE ANALOGY 
A. Time independent Hamiltonians: detectors 

Contrary to what could happen in a general relativis- 
tic context, in non-relativistic quantum mechanics, the 
notion of ground state (vacuum) and excited states (par- 
ticles) is globally defined (see, e.g. , [5S^). Moreover, if 
we consider static configurations {i.e. time- independent 
Hamiltonians), any energy eigenstate is preserved by the 
unitary time evolution and, thus, particle creation effects 
cannot take place. 

In order to overcome this obstacle and make radiation 
effects possible, it is necessary that the notion of particle 
that the detector measures does not correspond to the 
excitations of the Hamiltonian. This is not so unnatural 
since detectors are local objects, and therefore, there is 
no reason to think that they are be able to capture the 
precise structure of the eigenstates of the global Hamil- 
tonian. 

In particular, we need to introduce a notion of particle 
(with an associated momentum k and energy e(fc)) given 
in terms of some ladder operators rjk. Once the notion 
of particle that the detector measures is established, we 
can determine the momentum distribution (number of 
particles with momentum k) of the ground state of the 
system 

nik) = {GS\4rj,\GS), (1) 

and see whether it follows a thermal distribution. 

The notion of particle rjk, together with its associated 
momentum k and energy e(fc), define a test Hamiltonian 
HjjDW = X]fc ^i^)''lknk- The momentum distribution ([T]) 
that the detector measures is given by the overlap be- 
tween the ground state of the system and the eigenstates 
states of the test Hamiltonian Huj^y/- 



^ It is not obvious that it is impossible to reconcile these models 
with a gauge fixed version of Lorentz and difTeomorphism in- 
variant models. It is likely that in the continuum limit, these 
would be described by effective theories of the class of theories 
recently introduced by Hoi'ava |18| and extensively discussed in 
the literature, as it happens in Causal Dynamical Triangulations 



B. Time dependent Hamiltonians: quenches 

Another possibility for simulating particle creation 
phenomena is to consider time dependent Hamiltoni- 
ans. If some parameter of the Hamiltonian changes fast 
enough, the adiabatic approximation breaks down and 
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the system leaves its ground state. The populations of 
the excited states become non trivial and hence radiation 
is produced. 

An example of such a process is a time dependent tun- 
nelling amplitude between the sites of the lattice. This 
can be interpreted as a time dependent scale factor of the 
spatial geometry, giving the possibility to study, along 
the very same lines described for a black hole configura- 
tion, analogue cosmological particle production phenom- 
ena. 

A particularly interesting case of a time dependent 
Hamiltonian is the quench setting. A quench experiment 
consists of preparing a quantum system in the ground 
state of a certain Hamiltonian, and, suddenly, changing 
the Hamiltonian such that the quantum state does not 
correspond to an eigenstate of the new Hamiltonian any- 
more. The system, then, is out of equilibrium an evolves 
non-trivially in time. The goal of quench experiments is 
to study the time evolution of systems and in particular, 
its possible relaxation to equilibrium. 

Let us notice that the process of equilibration after a 
quantum quench for the whole system is, strictly speak- 
ing, never possible. This is due to the fact that the sys- 
tem evolves according to a unitary evolution. Neverthe- 
less, what we mean by equilibration is the relaxation of 
the expectation values of a certain family of observables 
to constant quantities. The question whether the equilib- 
rium state after a quench corresponds or not to a thermal 
state has been studied in several works [?7H5T] . 

If the Hamiltonian is quenched to the Hamiltonian 
that defines the notion of particle in the detector setting 
HudWj then the momentum distribution of the system at 
any time after the quench is identical to the momentum 
distribution measured by the detector in the time inde- 
pendent case, and hence, both settings become equiva- 
lent. This is due to the fact that the observables ?7^?yfe 
commutes with the Hamiltonians considered before and 
after the quench. 

The equivalence between the two settings suggests the 
use of quench experiments to simulate particle produc- 
tion also in the detector framework. Furthermore, the 
fact that the momentum distribution does not change 
along the evolution dictated by Hjjdw allows to mea- 
sure it at the most convenient time. 



C. The model 

The Hamiltonian we will consider in the present paper 
is the following Fermi-Hubbard model: 



background graph F. The matrices A\j and Bi^ are, 
respectively, the adjacency matrix of F and its antisym- 
metrized form"*. The sum runs over all the Ny vertices 
of the graph F. The coupling J is the tunneling of the 
particles between two connected sites and A controls the 
strength of the Hamiltonian terms that do not conserve 
the number of particles. 

In quantum graphity, the role of curvature in a con- 
tinuous space-manifold is played by the connectivity of 
a dynamic graph, while the backreaction of matter on 
geometry, in graphity models, is controlled usually by a 
term in the Hamiltonian that annihilates particles and 
creates links of the graph^, with the idea that the larger 
is the number of particles in a region, the larger is the 
effect on the connectivity of the curvature. 

We derived this model as an effective description of 
Quantum Graphity in the limit of a very dilute mat- 
ter content and a weak backreaction term between the 
particles and the graph (see Sec. IV). This violation of 
the number of particles related to the connectivity of the 
graph is produced in Hamiltonian ([2| by the second term 
proportional to A. In a more rigorous way, Hamiltonian 
([2]) is exactly obtained from the quantum graphity Hamil- 
tonian after freezing the evolution of the graph and as- 
suming that this is in a superpositon state (see Sec. IV 
for details). 

Let us mention that a similar effective model has been 
studied both theoretically and experimentally in the con- 
text of optical lattices [21] in order to give experimental 
evidence to the Lieb-Robinson bounds. We thus expect 
that the results of our paper could be potentially tested 
experimentally, given that the experimenter should be 
able to measure our observables. 



D. Detector and notion of particle 

In order to study particle creation effects, we need to 
specify the notion of particle (and its corresponding en- 
ergy) that our detector will measure. In a continuous flat 
space, we can define particles as plane waves with a well 
defined momentum k and energy h^k"^ / (2m). We have to 
generalize this idea of particle to graphs. 

The simplest graph one can consider is the discretiza- 
tion of the line or a ring lattice (in which we have under 
control irrelevant IR divergences). In the ring graph in 
particle, it is natural to define the annihilation operator 
of a particle with momentum k as the discrete Fourier 



ff(F, A) = - J E 4?«la. + - <^ +h.c.) , 

(2) 

where ai^a\ are annihilation/creation fermionic operators 
that annihilate/create a particle in the vertex i of the 



* The antisymmetrization is arbitrary up to a certain extent. In 
fact, Bij = j)Aij , with <T{i,j) being any function such that 
""(^li) = """(ii*) s-iid |cr(i,j)| = 1 Vi,j. However this does not 
change the properties of the system as long the graph is chosen 
properly. 

In order that the Hamiltonian is Hermitian, the opposite process 
is also required. 
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transform of the annihiliation operators in position a,- , 



L-l 



3=0 



with L the number of sites of the ring. Notice that these 
are precisely the eigenmodes of the hopping Hamiltonian 



L-l 



L-l 



Hmng = - ^(ajflj+i + h. c) = ^ 4k)vl 



fc=0 



where e(fc) — — 2 cos(27rfc/L) gives us the energy of the 
mode k. In the continuum hmit, both the wave function 
and the dispersion relation are recovered in their stan- 
dard form. 

Thus, when other lattice configurations are considered 
(torus, cylinder, etc.), it is natural to define particles as 
the eigenmodes of the hopping Hamiltonian supported 
on a regular graph. This hopping Hamiltonian is written 
as 



UDW 



fc=i 



eikWkVk, (3) 



where A^^ ° is the adjacency matrix of the regular graph 

Fq. The eigenmodes of Hjjdw, labelled by k, and with 
energy e(fc), define our notion of particle. These are cre- 
ated and annihilated by the operators rjj, and rjk , and are 
the excitation the detector measures. 

Let us mention that while Tq is a regular graph, in 
general, the graph T that defines Hamiltonian (2]) does 
not have to, since it will have regions with a higher or a 
lower connectivity. Roughly speaking, we can construct 
r by taking the flat graph Tq, and adding or removing 
links to/from it. 



E. Calculation 

The aim of our calculation is to probe how many par- 
ticles (as the ones defined as in the previous section) will 
be present in the ground state of the Hamiltonian of ^ . 
Therefore, we need to compute 



n(fc) = {GS\iilvk\GS) 



(4) 



and check whether the ground state, with respect to this 
notion of particle, appears as a thermal state, i.e. the 
momentum distribution Q follows a Fermi-Dirac distri- 
bution. 

In order for the computation of Q to make sense, the 
Hilbert space of the full Hamiltonian and the one of the 
detector have to have some overlap (in terms of unitary 
mappings from one to the other) . It is sufficient that the 
two different Hamiltonians are defined as operators on 
the same Hilbert space, i.e. the graphs have the same 
number of nodes. 



The Hamiltonian ([2]) is a standard quadratic model, 
hence, it can be diagonalized as 



L-l 



(5) 



by means of a Bogoliubov transformation of the funda- 
mental particle operators, a|,aj. In turn, these are re- 
lated by another Bogoliubov transformations to the op- 
erators rj, . Then, the operators r;, r]^ will be connected 
to the ■(/': V'^ by the Bogoliubov transformation that is 
the composition of the Bogoliubov transformations that 
relate "0, t/i^ to a, a) and a, to rj, rj^ . It can be written 
formally as 



L-l 



rik 



(6) 



where Ukj and jS^j are the Bogoliubov coefficents. Parti- 
cle creation effects will be clearly related to them. 

Notice that in general the ground state of these 
fermionic models is not just the Fock vacuum of the 
quasiparticles (the eigenmodes of the given Hamiltonian) , 
but it is a Fermi sphere due to the appearance of some 
negative energy modes in the spectrum (actually, this is 
the reason why we are using a Fermionic model, since 
a Bosonic one would lead to instabilities). This implies 
that the momentum distribution n{k) will receive a con- 
tribution just by the very presence of real quasiparticles 
in the Fermi sphere. 



III. ANALYSIS 

In this section we present numerical and analytical re- 
sults which support the claim of having particle produc- 
tion, at least in the sense of detectors. The key ingredient 
is the mismatch between the eigenspaces of the Hamilto- 
nian governing time evolution and the eigenspaces of the 
Hamiltonian with which we set up the detector. 



A. Numerical results 

We study here three different lattice geometries: a 
cylinder, a torus, and a ring. We evaluate the popu- 
lations {rjljjk) in the ground state of the eigenmodes of 
the Hamiltonian Hjjjjw that defines our notion of par- 
ticles, that is, the Hamiltonian with A = 0. Finally, we 
check if the distribution of the occupations with respect 
to the energy is thermal. 

In Fig.jl] the populations obtained in the three geome- 
tries are plotted with respect to the energy for A = 0.3. 
We observe that the observable is close to thermal in all 
cases. In particular, both for the cylinder and the torus, 
the energy modes as measured by a Unruh-de Witt de- 
tector are populated according to a Maxwell-Boltzmann 



5 



(a) Cigar 

Log[n(k)] 




(c) Ring 

LogLl/n(k)-l] 




e(k) 



FIG. 1: Functions of the number of quasiparticles n{k) = 
ivivk), evaluated numerically, in the ground state of the sys- 
tem with A = 0.3 for three different lattices: (a) cylinder 
[20x20], (b) torus [16x16] and (c) ring [300]. 

distribution, 

nMi3(fc)oce-^^W, (7) 

and therefore in the plan (log n{k), e{k)) we see a straight 
line with slope —(3. For the case of a a circle, we observe 
a Fermi-Dirac distribution, 

""FDik) = ^^J^ ^ ^ , (8) 

and the logarithm of the populations log(l/n(A;) — 1) 
scales linearly with e{k). 

For different values of the parameter A and size of the 
system, we numerically evaluated the Bogoliubov trans- 
formation between the modes and their populations n{k). 
In all cases, there are striking numerical evidences that: 

• The temperature is well defined in the large lattice 
size limit (thermodynamic limit). 



• The temperature (inverse of the slope in the plots 
of Fig. [ij scales proportionally to A, i.e. 

T cx A. 

• The higher the A, the better is the alignment of the 
points on a straight line, that is, the more thermal 
the behavior of the radiation. 

In order to better understand the thermal nature of 
these distributions, we solve in the next section the 1- 
dimensional chain with periodic boundary conditions an- 
alytically. 

B. Analytical solution of the 1-dimensional lattice 

In the one dimensional case, the Hamiltonian ([2| be- 
comes 

L-l L-l 

H = - J ^ (al^am+i + h. c.)-f A ^ {ajnam+i + H. c.) , 

m— m— 

(9) 

where L is the length of lattice (chosen to be odd) and the 
periodic boundary conditions are imposed by the identi- 
fication of ao with ql- Let us point out that this Hamil- 
tonian is equivalent to the XY model in an external mag- 
netic field [5^ . 

In order to diagonalize this model, let us exploit its 
translational invariance and introduce the Fourier modes 

1 

^'^ = 7t E e-^'^'-a™ . (10) 

^ m— 

We note that according to this definition the mode rj-k — 
i]L-k- From now on, we will use both [—{L — l)/2, (L — 
l)/2] and [0,i — 1] as the range of k depending on the 
simplicity of the notation. In terms of these modes, the 
Hamiltonian ^ reads 




L-l 

+ A ^ (e-'^'^TjkilL-k + h.c^ . 

fc=0 

Reshuffling some terms and using the canonical anti- 
commutation relations, the Hamiltonian can be written 

as 

L-l 

H=Y, Hk)vim + i(j{k)if]kVL-k - vl-kVl) , (11) 
where 

e{k) = - Jcos(^^k^ (12) 
cr(fc) = -A sin (^fc) • (13) 
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Let us note that while e(fc) is a symmetric function e(fc) = 
e{L — k), a{k) is anti-symmetric, cr(fc) = —a{L — k). 

Thus, the Hamiltonian Q is decomposed in {L — l)/2 
independent Hamiltonians 

L-l 

H=Y,Hk, (14) 

fc=0 

with 

Hk = e{k){r]lT]k+V-kV-k)+icr{k){r]lT]lk-V-kVk) , (15) 
or, in matrix notation, 

H.^(.-I .--.)(_f,i,lt',)(,Lj- <^'« 

1. Diagonalization 
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FIG. 2: Dispersion relation e(k) for different values of J and 
A; A = (continuous), |Jj > |A| (dashed), |J| = |A| (square 
function), |A| > |J| (points). 



First of all, let us note that 

e(fc) -\a{k) 
iCT(fc) -e(fc) 



(j{k)ay + e{k)(7^, (17) 



where the standard Pauli matrices. In the vec- 

tor space formed by the (1, cr^;, Cy, Cz), we have to rotate 
the vector (0, 0, a, e) in order to get a vector of the form 
(0, 0,0,w). Hence, one immediately realizes that the di- 
agonalizing matrix is a rotation of an angle 



6{k) — arctan 



e{k) 



(18) 



around the X axis. 



U = exp ( i-fTj. 



cos(6l/2)I-|-isin(6'/2)cr:^ . (19) 



From Eqs. (161 and ( 19 ), each Hamiltonian Hk can be di- 
agonalized with a Bogo iubov transformation of the form: 



V'fc = cos I — 1 rik - ism I — 1 r/i^ . 
The Hamiltonian in terms of its eigenmodes reads 



with 



ujk = sign(efe)Ve^(fcyTa^(fc) . 
The eigenvalues have then multiplicity two 

uj{k) =^ uj{L — k) ~ uj{—k) , 



(20) 



(21) 



(22) 



(23) 



that is, the energy of a quasiparticle with momentum k is 
the same as the energy a quasiparticle with momentum 
-k. 



2. Ground state 

For simplicity, we fix J = 1 and consider A > in 
units of J by rescaling the Hamiltonian. The ground 
state of the system with A = is a Fermi sphere, and the 
first excited states are gapless modes {e{k) oc {k — kp)). 
In the limit of large L, this is a conductor. This is a 
direct consequence of the Goldstone theorem: the Fermi 
sphere ground state spontaneously breaks the global U{1) 
symmetry of the original Hamiltonian, i.e. 

H{a,a)) = iJ(e'"a,e-'"a^). 

However, as expected from the fact that the A term 
breaks the global U{1) symmetry of the standard action, 
the first excited states are separated from the Fermi sur- 
face by a gap of 2A. In the limit of L large, and large A, 
this is an insulator. We can write then the ground state 
as 



\GS)= n v4io) 



(24) 



\k\<kB 



where |0)^ is the Fock vacuum for the ipk operators, and 
kp ~ (L — l)/4 is the the radius of the Fermi sphere 
(a;(fc) < for |fc| < kp). 

From Eq. ( p4| , we can derive the structure of the 
ground state in terms of the operators 77^ and (up 
to a global phase) 

\GS)=X{{sin{e{k)/2)~\co^{e{k)/2)r,y_^) 

k=Q 

(L-l)/2 

X n (cos(0(fc)/2)+isin((?(fc)/2)r;lr;L,)|O),, 



k=kp+l 



(25) 



and the occupation of the mode k reads 



(a) 



n(fc) := {GS\4rj,\GS) 



sm^{e{k)/2) if A: > fcj 



This piecewise function can be written in a more compact 
form by using the sine and cosine half angle formulas, 



n(fc) 



1 



2w2(fc) + e(fc)|a;(fc)| 



(26) 



Let us note that, fixed A, the momentum distribution 
n{k) only depends on its energy, since a and uj can be 
expressed in terms of e via cr^ — X^{1 — e^) and = 
a^ + e\ 

This has a nice interpretation. Due to the presence 
of the Bogoliubov transformation, the Fermi sphere ap- 
pears to be slightly depopulated, with the particles being 
radiated into the higher energy levels. 



3. Analogue thermal particle production 

In order to check whether the momentum distribution 
of the radiation measured by our detector in Eq. ( 26 ) 



is thermal, we have to compare it with the Fermi-Dirac 
distribution defined in Eq. 

In Fig. [sj^a), the quantity log(ri(fc)~^ — 1) is plotted 
with respect to £{k) for several values of A. If the radi- 
ation is thermal, these points should be fit by a straight 
line that passes through the origin an has slope p. This 
is precisely the behavior observed. In Fig. [sjjb), the real 
population (continuous line) of each quasiparticle with 
momentum k is compared to the population provided by 
the Fermi-Dirac distribution (dashed-line). Both distri- 
butions are very close to each other. 

In order to get an analytical expression for the tem- 
perature of the radiation, let us notice that, if the distri- 
bution is a Fermi-Dirac, then 



dn 
d^ 



(27) 



e=0 



In our case, the n{k) distribution is very close to a Fermi- 
Dirac, and therefore, the inverse temperature given by 
Eq. (27) can be approximated from Eqs. (12 1 and (26), 



(28) 



matching the expectation that T cx A. 



C. Comment on Thermalization 

From the results of the last section we conclude that 
the momentum distribution of the particles that the de- 
tector would measure is very close to thermal. Let us 
note that, although the observables seem those of a ther- 
mal system, this does not mean that system is thermal. 



Log[l/n(k)-l] 




FIG. 3: (a) Linear fit of the points (log(n(fc)"^ - l),e{k)) 
where n(k) is the momentum distribution of the modes rjk in 
the ground state of the Hamiltonian ([9| and e(fc) is the energy 
of each mode. Its alignment implies that the occupations of 
the modes k follow a Fermi-Dirac distribution with inverse 
temperature corresponding to the slope, (b) Comparison be- 
tween the real (continuous line) and the Fermi-Dirac (dashed 
line) distributions. 



There are other observables (apart from the number of 
particles T]l,r]k considered previously) for which the sys- 
tem will not look like thermal. But then, why does the 
momentum distribution n(fc) follow a thermal distribu- 
tion? 

In order to understand the origin of this thermal be- 
havior, it is useful to write the state of the system in 
terms of the eigenbasis of the Hamiltonian. In Ref. |33| . 
thermalization of local observables is explained by as- 
suming a narrow energy distribution together with typi- 
cality arguments [31] . In Ref. [551 ISS] , thermalization is 
a consequence of a narrow energy distribution and the 
Eigenstate Thermalization Hypothesis. In both cases, 
the temperature of the equilibrium state is determined 
by the position of the microcanonical energy window in 
the spectrum of the Hamiltonian. Determining the en- 
ergy distribution of the state could explain the approxi- 
mated thermalization observed by means of one of these 
mechanisms. 

Let us write the Hamiltonian (|3| that defines our no- 
tion of particle in its spectral representation for the 1- 
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dimensional case case, 

L-l 



H, 



UDW 



E^W'^I^fe^ E ^m|m)(m|, (29) 

m6{0,l}^ 



k=0 



where its eigenbasis {|ni)} is the Fock basis of the k 
modes, 



L-l 



(30) 



fc=0 



with m tlic bit-string of occupations with components 
"^fc = {0; 1} V k. The energy E^^ of each eigenstate |m) 
is given by 



L-l 



(31) 



k=0 



The first sum in Eq. (29 1 runs over the L different k- 
modes, while the second one is performed over all the 2^ 
elements of the eigenbasis of the Hamiltonian. The state 
of the system in this basis becomes: 



\GS) = J2 ^-1"!)' 

mG{0,l}^ 



(32) 



with = {m\GS). These coefficents can then be de- 
rived from Eq. ([25|) and written as 



cos(6i(/c)/2)'"'=J„ 



= e''^'" ]^sin(6i(fc)/2)^ 

fe=i 



(L-l)/2 



Y[ sin(0(fc)/2)'"'=cos(0(fc)/2) 



fc=fci!.+l 



where tpm is a phase given by 

/ kF 



fc=l 




(33) 



(34) 



The Kronecker delta Smk,m^k comes from the fact that 
\GS) is a superposition of only those Fock states with 
both modes k and —k either simultaniously occupied or 
empty (see Eq. (25)). This implies that only 2^/^ coefh- 
cents among the 2^^ possible ones are non-trivial. 

In Fig. |4j the absolute value of the coefficients of 
the initial state in the Hamiltonian eigenbasis given by 
Eq. ( p3l ) are plotted with respect to the energy. We ob- 
serve that they decrease exponentially. The coefficient 
of the exponential decay can be numerically determined, 
and we have to check that for any value of A it coin- 
cides with /3/2. The initial state can be then very well 
approximated by 



\GS) ^ Yl e'^-e- 



I3/2E„ 



m 



(35) 



where the sum only runs over the set of bit-strings M 
defined by 



M ^ {me {0, 1}^ I TUk = m^k V k] 



(36) 



Any observable A which equilibrates, will equilibrate 
towards its time average expected value. 



A — lim — 



dt{A{t)) = Tr(ia;) 



(37) 



where lo is the time average state defined by 

uj^ \im - dte"'^^°«'-*|G5)(G'5|e'^^°»'*. (38) 



Notice that w is the apparent equilibrium state. We say 
"apparent" because the system is never in ui, however, 
all those observables which equilibrate, do it towards the 
expected value that would be measured if the system was 
in uj. Let us point out that in the case that the Hamilto- 
nian has no degeneracies, w is the completely dephased 
state in the Hamiltonian eigenbasis and it maximizes the 
von Neumann entropy given the constants of motion |37j . 
Thus, the time average state uj reads in our case 



E 

me A/ 



Zm 



|m)(m| , 



(39) 



meM ' 



with Zm = ^ 

Now we can understand the origin of the Fermi-Dirac 
distribution of the populations of each mode 



n{k) = TT{filr]kUj) 



E 



1 



me{M\mk = l} 

where we have used that 



Z 



M 



Let us remark that uj in Eq. ( 39 1 is not a thermal state. 



since the sum does not run over all the eigenstates of the 
Hamiltonian, but only over those that have ruk = m^k- 
Nevertheless, in practice, all the observables that do not 
contain correlations bewteen the modes k and —k will 
give the same value as if the system was in a Gibbs state. 

In conclusion, we have seen that the thermal spec- 
trum of the radiation is a consequence of the exponen- 
tial decay of the energy distribution, that is, of the fact 
that the coefficients of the initial state written in the 
Hamiltonian eigenbasis decay exponentially (see Fig. |4] 
and Eq. (35 1). This exponential decaying is not generic 



me A/ 



at all, in the sense that, most of Bogoliubov transforma- 
tions relating the eigenbasis of two different Hamiltonians 
(even if they have the same locality structure) will not 
produce it. This represents another mechanism towards 
thermalization, where no bath is necessary and where the 
temperature is not given by a macroscopic energy scale 
(the position of the microcanonical window of energies 
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Log[|Cn| 




En 



FIG. 4: Coefficients of tfie state of tlie system in tlie Hamil- 
tonian eigenbasis witli respect to tlie energy in a logscale in 
tlie vertical axis. Their alignment in a straight line shows an 
exponential decay with the energy. 



in the spectrum) but by the intensive parameter which 
is quenched. It is an open question beyond the scope 
of this article to understand what properties two Hamil- 
tonians have to share in order for the coefRcents of an 
eigenstate of one Hamiltonian, spanned in the eigenbasis 
of the other, to decay exponentially respect to the energy. 



IV. APPLICATION TO TRAPPING SURFACES 
IN QUANTUM GRAPHITY 

In this section we show how the model considered 
in the present paper can be derived from Quantum 
Graphity. Moreover, we will apply the procedure to the 
trapping surfaces considered previously in [23] . 



Derivation 



nians read 

-fflinks 



-UT.-h)^ (41) 

(id) 



Hbr = a ^ (^S-^ j-^ (E) ala] + S+ ..^ (g) a,aj) , (43) 
(hi) 

where Pij = is the projector onto the edge 

S± = |l)(0|jj and S'^ = |0)(l|jj, and U, t, and A set the 
speed of the oscillations of the links, the tunneling rate 
and the backreaction strength respectivelly. Notice that 
in the Hamiltonian for the matter -ffparticics the interac- 
tion among the particles has been neglected, but natu- 
rally it can be included. However, for the phenomena we 
are interested in, we do not really need interactions, and 
we can work with free fields. First of all, let us assume a 
situation in which matter degrees of freedom do not af- 
fect the graph, i.e. their coupling to the dynamical graph 
is small compared to hopping one. This means that the 
edges of the graph can be considered to be constant in 
time, and thus a fixed graph in the first approximation. 
However, we will still retain the quantum structure of 
the graph, i.e. links can be still in a superposition of on 
and off states. In this approximation, we can derive the 
Hamiltonian ^ from Quantum Graphity. 

In order to see this, we consider the matrix-elements 
of the Hamiltonian ( |40[ ) applied to a particular quantum 
graph, and realize that these are the same as those of the 
Hamiltonian in Eq. ([2]): 

((l),n\H\(f>n) = (((/),„| (g> (l/'spacoD-fftotal {\<Pn) IV'space)) , 

(44) 

with the quantum state for the links fixed to 



(45) 



One of the most problematic features of Quantum 
Graphity models is their complexity. The dynamics of 
the graph makes the model hard to solve analytically. 
Thus, approximate techniques to understand the behav- 
ior of the system are required. 

First of all, we briefly describe the structure of the 
Hubbard model on a dynamical graph introduced in 
Ref. [H]. The Hilbert space of the system is 'Hnnks ® 

n particles 

with "Hlinks — ^i^l " 'Hparticlcs — 

and Hi the local Hilbert space of one site. Its 
Hamiltonian can be written as 



total 



links 



particles 



BR 



(40) 



where iJunks describes the dynamics of the graph (space), 
iJparticics is the Hamiltouian for the particles (matter), 
and 77br corresponds to the backreaction interaction be- 
tween space and matter. More explicitly, these Hamilto- 



where Q is a Hadamard gate, Q\0) = 1/72(10) + |1)), 
and Aij is the adjacency matrix of the graph. 

The state iV'space) is then a product state of all the 
links in the state |0) except those links that we would 
like to have active in the effective model, which are in 
the superposition 



l+) = -^(|o) + |i)) = g|o) 



(46) 



In order to prove Eq. (44) it is only necessary to realize 
that 



(+|P|+) = (+|5±|+) = ^ 
(0|P|0) = (0|5±|0) = 0. 



1 



(47) 
(48) 

We argue then that the Hamiltonian of Eq. ^ can be 
considered as an effective model in the case in which the 
graph is in a equal superposition of on and off states. 
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Some comments are in order. First of all, let us no- 
tice that the bosonic model is unstable in the thermody- 
namic limit in which the number of nodes of the graph 
goes to infinity. This is due to the fact that the energy 
can be decreased by increasing the number of bosons in 
the graph. Morever, the model is reliable as an effec- 
tive model of Quantum Graphity as long as the num- 
ber of particles Np < 2L + Ni, where L is the number 
of links and Ni is the number of particles in the states 
iV'totai) = lipgraph) IV'part)- The rcason is that in Quan- 
tum Graphity the quantity Q — N — 2L is conserved. 



B. Unruh— de Witt detectors 

In order to apply the procedure we described in the 
first part of the paper to Quantum Graphity, it is neces- 
sary to introduce the proper interpretation to the Unruh- 
de Witt observables introduced previously. We use the 
fact that particle states are identified by the behavior 
of suitably defined detectors, i.e. devices that count the 
number of excitations associated to the measured quan- 
tum state. As in the case of QFT in curved spaces, these 
detectors will be associated to an observer, and hence an 
Hamiltonian with respect to which states are classified. 
In the case of QFT in curved spaces, particle production 
as the Unruh effect can be purely kinematical phenom- 
ena associated to the mismatch between the notion of 
particle associated to different observers, possessing, ef- 
fectively, different Hamiltonians^ . Here we consider the 
ideal experiment of an observer that constructs his detec- 
tor in a region where the background graph is essentially 
regular (flat), and hence moves towards a region of high 
connectivity, interpreted as a black hole analogue. The 
readings of these detectors will be then influenced by the 
change in the notion of particle, ultimately associated to 
the structure of the graph. The calculations, then, are 
the same as those in the first part of the paper, although 
the quench involves also the structure of the graph, apart 
from the parameter A. 



C. Difficulties with the particle production 
interpretation 

In [331 US] it has been shown that certain graphs con- 
figurations might work effectively as black holes. It has 
been shown that curvature is related to variations in the 
connectivity of the graph. A rather extreme case is when 
the otherwise random graph contains a highly connected 
region, i.e. a subset of nodes which are connected by 
a link with almost any other node of the subset itself. In 



The evolution of an uniformly accelerated observer is controlled 
by the boost Hamiltonian. 



a certain limit (related to the size of 3§) this region be- 
comes a trapped region, i. e. a region from which particles 
cannot escape, once they enter it. 

As discussed in [23], this phenomenon is rather pecu- 
liar with respect to the usual classical intuition in terms 
of bundles of geodesies. Indeed, it can be shown that 
the lower energy eigenstates correspond to the case in 
which the particles are trapped into this region. More- 
over, the gap between the ground state and the excited 
states (with zero angular momentum) is proportional to 
the local degree (connectivity) of the graph. The size of 
the gap is finite for finite graphs, and hence these regions 
are not completely trapping. However, the larger is the 
size of this region (and thus the gap), the larger energy 
needs to be transferred to a trapped particle to kick it 
out of it. 

We are interested in the following question: can we see 
Hawking radiation? Unless we use this quench formalism, 
this is rather unlikely. 

One way to see this is to say that, while the region ^ 
behaves as a trapping region, it is not associated to an 
ergoregion where negative energy modes can be stored. 
A more precise statement involves the unique well de- 
fined notion of ground state and particles so far available 
in these models. These are limitations typical of any or- 
dinary non-relativistic approach to this problem and as 
such, also Quantum Graphity suffers of it. 

However, we will ask the following practical question. 
How will a particle detector react when lowered near the 
region ^? An important point to make explicit here 
is that the detector will work as an external apparatus, 
i.e. is not assumed to be built in terms of the same mat- 
ter fields living on the graph. If this black hole configura- 
tion were realized in a lab, the detector we are referring 
to will be a scanner, a particle counter operated from 
the experimenter in the lab, who would just count the 
clicks corresponding to suitably defined notion of parti- 
cle. In turn a particle is defined as an excitation above 
the ground state of a given Hamiltonian. 

Formally, the setup described in the previous sections 
apply straightforwardly to the trapped surface lattice 
described above. For increasing values of A, the loga- 
rithm of the population number becomes more and more 
a straight line with negative slope. This result is consis- 
tent with the findings in the 1-dimensional model. We 
thus expect that this phenomenon might be, to a cer- 
tain extent, independent from the dimensionality of the 
lattice or from its connectivity. 

We performed simulations for the lattice with cylindri- 
cal topology and with a complete graph attached to one 
end of the cylinder. This graph has been as a model of 
trapped surface in [23|. We report the results here for 
simplicity, as we do not believe the plots have particular 
relevance, other than showing a particular slow conver- 
gence. The parameter M which is the number of nodes in 
the angular coordinate of the cylinder, is also the size of 
the complete graph. We observe that for increasing val- 
ues of M, the temperature obtained from the fit of the 
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logarithm of the population number is proportional to the 
value of A. There is a reason for this, which we believe is 
easier to understand in the effective 1-dimensional model 
obtained in In this case, the Hamiltonian for the 

modes with zero angular momentum is given, for A = 0: 

i?, 

iJeff = - ^ Ma\a,+i + h. c.) (49) 

where Jq ^ MJ and Jt ^ J for i > 1. When M > 1, 
the effective Hamiltonian of the model is the one with 
a node disconnected from the rest of the graph. In 
the original 2-dimensional model, this means that the 
complete graph can be considered as disconnected from 
the rest of the graph and the ground state becomes the 
product of the ground state for the complete region and 
the one for the cylindrical lattice, which then reduces 
the analysis to the one of the cylindrical lattice. Thus, 
the complete graph, as a matter of fact, does not change 
the temperature of observed radiation but, indeed, the 
more the region is trapping, the more the temperature 
converges to the value of the flat one. This is in fact a 
negative result, as it means that the radiation observed 
by the Unruh~de Witt detectors are not related to the 
trapping region, but only the quenching of the parameter 
A. 



V. CONCLUSIONS 

In this paper we investigated the possibility of using 
quantum quenches for simulating effects of quantum field 
theory in curved spacetime. This research was inspired 
by ideas arised in the context of Quantum Graphity mod- 
els, but can be placed also in the context of analogue 
models. 

A certain number of analogue models have been pro- 
posed in order to simulate gravitational phenomena or 
effects of quantum field theory in the presence of curved 
geometries, e.g. Hawking radiation and cosmological par- 
ticle production. In our case, our proposal is different 
from previous models since we do not rely on the ap- 
pearance of effective causal horizons, but rather on a dif- 
ferent origin of the mismatch between notion of particles 
as they are defined for different detectors. 

In fact, we have shown that by considering Hamiltoni- 
ans which do not conserve the particle number operator 
we can simulate approximate thermal particle production 
by means of a Quantum Quench. Quantum quenches 
are very interesting for many reasons. First of all, for 
many quantum mechanical systems it was shown that 
many observables thermalize in the long time regime af- 
ter a quantum quench. The mechanism is still poorly 
understood, even though it is thought that this might 
be related to the eigenstates thermalization hypothesis 
[3Sll3n]- The setting we proposed is related, as we have 
seen, to the quantum quench of one parameter in the 



Hamiltonian. We proposed to consider first the ground 
state of a Fermi-Hubbard model and then to perform 
a sudden quench of a parameter associated to a par- 
ticle number non-conserving term in the Hamiltonian. 
We have shown analytically and numerically that the 
expectation value of Unruh-de Witt detectors is very 
close to thermal. This can be achieved by relating the 
quenched and unquenched Hamiltonian through a Bo- 
goliubov transformation. 

Some comments are in order. First of all, the expec- 
tation values of the number of particle in each mode do 
not evolve in time, thus their value just after the quench 
is the same they would have long after it. This situation 
is similar to the one appearing in the context of quan- 
tum quenches, long after the quench, with the important 
difference that, contrarily to time-evolving observables, 
these can be measured right after the quench and do not 
need to wait the "relaxation" of the system. 

The search for thermality is a very common aspect in 
quench settings. While the total system is in a pure state, 
the expectation value of some observable are very close 
to thermal in the long time regime. This is confirmed by 
our calculations and numerical simulations, as the inspec- 
tion of the occupation number of modes corresponding 
to particles detected in analogue detectors. This could 
represent a new connection between the broad area of 
quantum systems in the laboratory and simulators for 
gravity. 

However, the direct application to the case of black 
holes, as we have argued, is not appropriate, given that 
the same particle creation effect is present in backgrounds 
that do not possess obvious causal boundaries. 

Indeed, when analyzing the case of black hole-like 
graphs appearing in Quantum Graphity, to see whether 
they could be emitters of Hawking radiation, we could 
show that the temperature of the Planck distribution did 
not depend on the size of the trapped-surface, but indeed 
converged to a value independent from the geometrical 
properties of the trapped surface. We gave an explana- 
tion for this, based on the fact that the effective coupling 
between the trapping region and the rest of the graph 
decreases with the size of the trapped surface. 

Nonetheless, along the lines of [3], the setup that we 
described in the present paper could be used to simulate 
cosmological particle production, making the parameters 
of the Hamiltonian functions of time, controlled by the 
experimenter. 
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